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Self-assembly dynamics in binary surfactant mixtures and structure changes of lipid vesicles in- 
duced by detergent solution are studied using coarse-grained molecular simulations. Disk-shaped 
micelles, the bicelles, are stabilized by detergents surrounding the rim of a bilayer disk of lipids. 
The self-assembled bicelles are considerably smaller than bicelles formed from vesicle rupture, and 
their size is determined by the concentrations of lipids and detergents and the interactions between 
the two species. The detergent-adsorption induces spontaneous curvature of the vesicle bilayer and 
results in vesicle division into two vesicles or vesicle rupture into worm-like micelles. The division 
occurs mainly via the inverse pathway of the modified stalk model. For large spontaneous curvature 
of the monolayers of the detergents, a pore is often opened, thereby leading to vesicle division or 
worm-like micelle formation. 



I. INTRODUCTION 

Amphophilic molecules such as lipids and deter- 
gents self-assemble into various structures in aqueous 
solutions 1 - - — . A single type of simple surfactant can 
form spherical and cylindrical (worm-like) micelles, bi- 
layer membranes, inverted hexagonal structures, and in- 
verted micelles depending on the relative size of their 
hydrophilic parts 1 -^. Aggregates of different types of 
surfactants can exhibit more complicated self-assembled 
structures such as undulated cylinders, octopus-like mi- 
celles, and vesicles connected with worm-like micelles 9 -*^. 
In mixtures of lipids and cone-shaped surfactants with 
a relatively large hydrophilic head, the edges of bilayer 
membranes are stabilized by surfactants, and a bilayer 
disk called bicelle is formed^—. The bicelle is an excel- 
lent medium to investigate protein functions, since mem- 
brane proteins can maintain their functionality in the 
bicelles&ii. Bicelles can align in magnetic fields and are 
used for NMR studies 1 ^. The structures of several pro- 
teins have been determined by the bicelle method&i 3 -. 

The dynamics of micelle growth and subsequent vesi- 
cle formation have been studied experimentally by 
time-resolved scatting techniques (light, neutron, and 
X-ray)^~— . Mixtures of lipids (lecithin) and natu- 
ral detergents (bile salts) 14 ' 17 or cationic and anionic 
surfactant o 15 ' 18 are often used in experiments, since 
the size and shape of the surfactant aggregates can be 
changed by the concentration variation. Bile salts are 
also typical detergents that are used to produce bi- 
celles. The self-assembly of lipids into vesicles in single- 
component lipid systems has been studied using simu- 
lations of a meshless membrane model with and with- 
out hydrodynamics 1 ^. However, in contrast, the self- 
assembly dynamics of lipids and detergents to form bi- 
celles has not been completely understood. The first aim 
of this study is to investigate the self-assembly dynamics 
of binary mixtures of lipids and detergents. 

When a detergent solution is added to a liposome sus- 
pension, the liposomes can dissolve. Experimental stud- 
ies have revealed that liposomes take various types of 



solubilization pathways depending on the types of lipids 
and detergents 2 ^ - — ; these including rhythmic shrink- 
age, bursting, budding, fission, peeling, and inside-out 
inversion. In particular, vesicle fission (or division) is 
considered an important dynamics since vesicle growth 
and division are two fundamental processes of the self- 
reproduction. Recently, the self-reproduction of vesi- 
cles has received growing attention in terms of both 
experiments 2 ^— and theor y 31 i 32 , since it is one of the 
key functions of protocells. To our knowledge, surfactant- 
adsorption-induced vesicle division has been numerically 
investigated only in a few simulations carried out by 
Markvoort and coworkers^ 3 .. The second aim of this 
study is to systematically investigate the dynamics and 
mechanism of detergent-adsorption-induced vesicle divi- 
sion. 

Membrane fusion and fission are key events in vari- 
ous intra- and intercellular processes in living cells, and 
they have been intensively investigated 3 ^—. The fis- 
sion dynamics of a bilayer membrane has been sim- 
ulated for single-component 3 - 9 - - — and two-component 
membrane s 33 ' 42 ' 43 . Membrane fission occurs via path- 
ways opposite to those of membrane fusion via a stalk 
intermediate. In one of the fission pathways, the for- 
mation of a trans-monolayer contact 4 ^ induces the stalk 
formatio n 33 ' 40 , and in the other, stalk formation is in- 
duced by opening of a small por o 39 ' 42 . The control pa- 
rameters that determine these fission and fusion path- 
ways have not been well understood so far. 

In order to simulate bilayer membranes on longer and 
larger scales, various types of coarse-grained molecular 
models have been proposed (see review articles 4 ^— ). 
In the present study, we employ one of the solvent-free 
molecular models^. In this model, membrane proper- 
ties such as the bending rigidity and spontaneous cur- 
vature of the monolayer can be varied over wide ranges. 
In our previous study-^, we found that the rupture of 
two-component vesicles can lead to the formation of var- 
iously shaped micelles and vesicles including bicelles and 
octopus-like micelles. We apply the same model to the 
molecular self-assembly into bicelles and detergent ad- 
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sorption to vesicles in this paper. 

In Sec. HH the surfactant model and simulation method 
are briefly described. In Sec. Mil self-assembly from 
isolated molecules into bicelles and vesicles are investi- 
gated. In Sec. HV] detergent-adsorption-induced topolog- 
ical changes, vesicle division and rupture into worm-like 
micelles are described. The summary is provided in Sec. 



II. MODEL AND METHOD 

We employ a solvent-free molecular mode l 50 ' 51 to simu- 
late two-component surfactant mixtures. One type (type 
A) of molecules is considered as the detergent molecule, 
being conically shaped with a relatively large head (the 
spontaneous curvature of the monolayer Co > 0). The 
other type (type B) is considered as a lipid with a cylin- 
drical shape (Co = 0). The surfactant model and simu- 
lation method are only briefly described here, since the 
details of the model can be found in Refs. [50ll5ll 

Each (i-th) molecule of the types A and B has a spher- 
ical particle with an orientation vector u^, which rep- 
resents the direction from the hydrophobic to the hy- 
drophilic part. There are two points of interaction in 
the molecule: the center of a sphere rf and a hydrophilic 
point r c t = r| + u^er. Type A molecules are considered 
for < i < N A and type B for N A < i < N = N A + N B . 
The molecules interact with each other via the potential 
given by 
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denotes the thermal energy. 

An excluded volume of the molecules with a diameter 
o and an attractive interaction between the molecules are 
given by the first and second terms in Eq. ([T]), respec- 
tively. The two types of molecules can have different at- 
tractive strengths: e% = £ att for < i < Na and £j = £ att 
for Na < i < N. The potential U a tt(pi) is given by 

U & M = 0.25 ln[l + exp{-4(^ - p*)}} - C, (2) 

with C = 0.25 ln{l + cxp(4p*)}. The local particle den- 
sity P i is approximately the number of particles in the 
sphere with radius r a tt, 
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where f C ut(r) is a C°° cutoff function^, 

/-«={7 {Ao(1+w ^ )} ^>;i;;J w 

with n = 6, A = ln(2){(r cut /r att ) n - 1}, r att = 1.9a, 
and the cutoff radius r cu t = 2.4o\ 

The third and fourth terms in Eq. (JTJ) are discretized 
versions of the tilt and bending potentials of the tilt 
mode l 53 ' 54 , respectively. A smoothly truncated Gaussian 
function^ is employed as the weight function 

Wcv (r) = { ^ ('/tT'-i ) ( r < r ") (5) 
[0 (r > r cc ) 

with n = 4, r ga = 1.5(7, and r cc = 3c. The spontaneous 
curvatures of the monolayer membranes are varied by 
parameters C b A d and C b B d : C^ = (C bd + C» d )/2, C b 



bd 



C^ for < i < N A , and C bd = C£ d for N A < i < N. 

The last term in Eq. ([T]) represents the repulsion be- 
tween the different types of molecules as a monotonic 
decreasing function: Uab{t) = Aif cut (r) with n = 1, 
-4o = 1, r-cut = 2.4cr, and A x = exp[tj/(r cut - a)}. 

The NVT ensemble (constant number of molecules 
A, volume V, and temperature T) is used with periodic 
boundary conditions in a cubic box with side length L. 
Undcrdamped Langevin equations [Brownian dynamics 
(BD)] are used for time evolution of the molecules. The 
motion of the center of the mass rp = (rf + rf)/2 and 
the orientation are given by 
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where m and / denote the mass and the moment of in- 
ertia of the molecule, respectively. The forces are given 
by fi G = —dU/drf and f/ = —dU/dui with the perpen- 
dicular component = a — (a ■ Uj)iii, and a Lagrange 
multiplier Al is used to ensure that uf = 1. The fric- 
tion coefficients and £ r and the Gaussian white noises 
g~ (i) and g£(t) obey the fluctuation-dissipation theorem. 

In this paper, we use fcbend = fctilt = 8, p* = 14, 
e att = 2, and C bd = 0. The type B (lipid) molecules 
form a fluid bilayer membrane with n/k^T = 44, and 
they have very low values of CMC. No molecules are 
detached from a single-component membrane of type 
B in typical simulation time scales. For the type A 
(detergent) molecules, the attraction range is set to 
0.6 < e A tt < 1.5 so that these molecules have a finite 
CMC range ~ 10~ 4 /cr 3 to 10~ 2 /cr 3 . The spontaneous 
curvature Co can be varied by the parameter Cbd as 

C 21 {fcbend/(fcbend + httt)}C h d/<j£^- At fc b end = &tilt, 

Co — Cbd/2cr. For type A molecules, the curvature pa- 
rameter range is set to C A d = to 1.2 (Cq — to 0.6er). 
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FIG. 1: Dependence of (a) the lateral diffusion coefficient D 
and (b) the effective line tension F of the membrane edges on 
£ab at TVa = A^b- (a) The symbols with lines represent data 
for C^d = (o), 0.4 (□), 0.8 (o), and f .2 (A) at e^t = 2. The 
diffusion coefficient D is normalized using the diffusion coef- 
ficient Do = c 2 / T ° °f an isolated molecule, (b) The symbols 
with solid lines represent data at e^tt — 2. The symbols with 
dashed lines represent data at e^ tt = 1.5 (V) and 1 (x). 



ccntrations of the type A and B molecules are given by 
p A = 0.0156/ct 3 and p B = 0.0039/ct 3 , respectively. Be- 
fore the execution of the product runs, the molecules are 
equilibriated at £ab = 10 so that no type A molecules are 
absorbed into the vesicle. At the beginning of the prod- 
uct runs (t = 0), £ab is changed to a negative value. The 
negative values of £ab are used to generate attraction be- 
tween type A and type B molecules. Positively-charged 
surfactants and negatively-charged surfactants can have 
similar attraction, and this attraction enhances absorp- 
tion of the type A molecules into the vesicle. 

The results are displayed with a length unit of er, an 
energy unit of k B T, and a time unit of t = ^a 2 /UbT. 
The error bars of the data are estimated from the stan- 
dard errors of ten independent runs. The time unit is 
estimated as tq ~ 0.1/xs from a ~ 2nm and the lateral 
diffusion coefficient ~ 10 _8 cm 2 /s for phospholipids^. 

Figure [1] shows the £ab dependence of the lateral dif- 
fusion coefficient D and the effective line tension V of 
the membrane edges. The quantities D and T are esti- 
mated from the mean square displacement and pressure 
tensor for Na = Nb, respectively^. As £ab decreases, D 
decreases and subsequently exhibits a jump into zero at 
^ab ~ — 6. At £ab % — 6, the membrane is in gel phase. 
Here, we mainly use the parameter region of fluid mem- 
branes. The line tension V of the edges of two-component 
membranes depends on £ab, since the positive sponta- 
neous curvature of type A molecules can reduce the free 
energy of the edges. As £ab increases, V decreases and 
reaches a minimum value at the phase separation point. 
For completely phase-separated membranes, T increases 
as the sum of the two line tensions: T = Fa + Tab^ [see 
the blue line with squares in Fig. El))]. With decreas- 
ing £^j. t at C^ d = 0.8, the line tension T decreases, and 
subsequently, the bilayer edges become unstable, thereby 
leading to the formation of worm-like micelles. 



A single-component vesicle is ruptured into worm-like mi- 
celles when C > 0.8 at e att = 

In order to investigate molecular self-assembly (Sec. 
IIIip . wc set randomly distributed gas states at £* tt = 
e att = 0-5 and C£ d — as initial states with TV = Na + 
Nb = 4, 000 and L — 40cr (the total concentration p = 
N/V = 0.0625/er 3 ). Subsequently, the simulations are 
begun with e£. t = 1.5 and C£ d = 0.8 (C£ = 0.4a). 

The repulsion between type A and type B molecules is 
varied from £ab = to 3. In tensionless membranes, the 
type A and B molecules are phase-separated at £ab ^ 1, 
and the line tension Tab between the type A and type B 
domains increases with increasing £ab^ Tab = 1-7 and 
3.5 at £ab = 1-5 and 2, respectively. At £ab ^ 2, flat 
membranes become unstable; the membrane is ruptured 
from the domain boundary, and subsequently, type A 
molecules form worm-like micelles. 

In order to investigate surfactant absorption into vesi- 
cles (Sec. IIV[) . the type A molecules are distributed 
outside of a vesicle composed of type B molecules for 
Na = 8,000, Nb = 2,000, and L = 80a. The con- 



III. SELF-ASSEMBLY INTO BICELLES 

Figures [2HS1 show the self-assembly dynamics of bi- 
nary mixtures of detergents (type A molecules) and 
lipids (type B). First, isolated molecules aggregate into 
small spherical micelles. These micelles fuse and grow 
into disk-like micelles. In a single-component system 
(N = Nb), large disk-like micelles of size i c \ > 1000 
close into vesicles via a bowl shape. In two-component 
systems, the type A molecules surround the rim of the 
bilayer disks of the type B molecules, and they can sta- 
bilize the disk-shaped micelles, bicelles (see Fig. [5]). As 
the number fraction N B /N of the type B molecules or 
£ab increases, larger bicelles arc formed. 

Figure [5] shows the time development of the mean clus- 
ter size (n c i) defined as 

Eoo -2 

n c i = ^oc : , (9) 
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FIG. 2: Sequential snapshots of self-assembly into bicelles at 
N A = N B = 2,000 and £ab = 1. (a-d) All molecules in the 
simulation are shown at (a) t = 0, (b) 1, OOOro, (c) 100, OOOro, 
and (d) 1,000, OOOro. (e-g) The largest cluster is shown at 
(e) t = 1, OOOro, (f) 10,000r , and (g) l,000,000r o . The red 
(yellow) and light blue (light yellow) hemispheres represent 
the hydrophilic (hydrophobic) parts of type A and type B 
molecules, respectively. 




FIG. 3: Sequential snapshots of self-assembly into a vesicle 
and micelles at N A = N B = 2, 000 and e A B = 3. (a-d) All 
molecules in the simulation are shown at (a) t — 5, OOOro, 
(b) 100, OOOro, (c) 803, 000r , and (d) l,000,000r . (e-g) The 
largest cluster is shown at (e) t = 5, OOOro, (f) 30, OOOro, and 
(g) 100, OOOro. 




where rii is the number of clusters of size i c \, and the nota- 
tion (...) denotes the average for ten independent simula- 
tion runs. We consider that molecules belong to a cluster 
when their distance is less than r at t = 1.9a. For binary 
mixtures, the growth of the clusters saturates as the disk 
rims are covered by the type A molecules [see Fig. EJb)]. 
In the single-component system, (n c i) increases linearly 
with time t until vesicle formation starts at t ~ 700, 000to 
[see Fig. [SJa)] . This behavior agrees with our previous re- 
sults obtained by simulations of the mcshlcss membrane 
model and by Smoluchowski rate equations^. In our pre- 
vious study, it was found that hydrodynamic interactions 
increase the rate of aggregation but do not change the 
qualitative dynamics of the self-assembly. Since hydro- 
dynamic interactions do not change the essence of the 
dynamics, in this study, we neglect them and use only 
BD. 

As £ab increases, the cluster size increases [see Fig. 
[S{c)]. This increase is caused by the formation of larger 
bicelles for £ab < 1 [compare Figs. [2] and B|a)], but 
by separation of the type A domain for eab > 2. At 
£ab = 2, the circumferences of the type B disks are 
not completely surrounded by type A molecules [see Fig. 




FIG. 4: Snapshots of self-assembled micelles for (a) eab = 
and (b) e AB = 2 at N A = N B = 2, 000 and t = 1, 000, 000r . 



HJb)], so that it is in a partial wetting condition of the 
type A domains at open edges of the type B domain s 51 ' 58 . 
At £ab = 3, the small micelles initially consist of both 
surfactants; however, the two surfactants eventually form 
completely separated aggregates. The type A molecules 
form spherical micelles and the type B molecules form 
bilayer disks or vesicles (see Fig. 
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FIG. 6: Time development of the mean asphericity {a sp } of 
clusters at Nb/N = 0.5, (a) Averaged for all clusters with 
n c i > 10 at £ab = 0, 1, 2, and 3. (b) Averaged for clusters 
with 10 < n c i < 50, 50 < n c \ < 100, and n c i > 100 at eab = 3. 
The error-bars are shown at several data points. 



FIG. 5: Time development of the mean cluster size (n c \). (a) 
Nb/N = 1. (b) Nb/N = 0.25, 0.5, 0.75, 0.9, and 1 at sab = 1. 
(c) £ab = 0, 1, 2, and 3 at Nb/N = 0.5. The error-bars are 
shown at several data points. 



The shapes of micelles are investigated by calculating 
a shape parameter, asphericity a sp , which is expressed 
as^ 

(A 1 -A 2 ) 2 + (A 2 -A 3 ) 2 + (A 3 -A 1 ) 2 
asp = 2(A 1+ A 2 + A 3 ) 2 ' (10) 

where Ai < A 2 < A3 are the eigenvalues of the gyra- 
tion tensor of each micelle. The asphericity is the de- 



gree of deviation from a spherical shape: a 



for 

spheres, a sp = 1 for thin rods, and a sp — 0.25 for thin 
disks. The time development of the mean asphericity 
(a sp ) = (J2 icl a spicini)/(J2 icl iciTk) is shown in Fig. |6] 
For £ab < 1) the increase in (a sp ) simply reflects an in- 
crease in the cluster size, since the a sp value of each size 
is almost time- independent: a sp = 0.1 (quasi-spherical 
shapes) and 0.22 (disks) for 10 < n c i < 50 and n c \ > 50, 
respectively. For £ab > 2, the cluster shape of each 
size changes in time [see Fig. EKb)]. Small micelles in 
the early stage, t < 100, OOOro, exhibit elongated shapes 
(a sp > 0.5), where one or two spherical type A domains 
are connected to a spherical type B domain [sec Figs. 



[31(e) and (f)]. Subsequently, the cluster growth results in 
the detachment of the type A domains, as the type A do- 
mains become sufficiently large to form stable separated 
micelles. This detachment yields a dip in the time evolu- 
tion curve of the mean cluster size (n c i) at t ~ 100, OOOro. 
After the detachment, the type B bilayer disks continue 
to grow, subsequently resulting in the vesicle formation 
[see Figs. [3jc) and (d)]. 

Next, we estimate the bicelle size in a monodisperse 
distribution. We assume that a bicelle consists of a bi- 
layer disk of the type B domain with radius i?dis — ho 
and a surrounding ring of the type A domain with ra- 
dius i?dis- The number of type A and type B molecules 
in the bicelle are expressed as nb c ,A = 47ri?dis/'o and 
"bc,B = 27r(i?dis ^ ^o) 2 / a o, respectively, where Iq and ao 
are the circumference length and area per molecule. The 
factor 2 is multiplied to the number of molecules because 
the membrane has two lipid layers. In the monodisperse 
distribution, the number of molecules are given by 



"-be, A = 47T, -j 




Qq N b 2h \ a N B 
ll N A l J ll N A 



8tt- 



ra Q i\ B n 1//10VI , 

r W^A + 7 --4l7 -) J ' ^Vb" 

and n bc ,B = n hc , A N B/N A for uq/Iq ~ 1 and h /l 



6 




FIG. 7: Sequential snapshots of vesicle division due to sur- 
factant adsorption at C^d — 0.8 and eab = —5. (a) t = 
(b) t = 8,550r . (c) t = 8,563r . (d) t = 8,570r . (a) All 
molecules in the simulation are shown, (b)-(d) Upper and 
lower snapshots show whole and cross-sectional-views of vesi- 
cles, respectively. 



1. For N B /N = 0.25, 0.5, 0.75, and 0.9 with a = a 2 
and Iq = h = a, the total number of molecules in the 
bicelle are 40, 90, 400, and 2,500, respectively. These 
estimations agree with the simulation results. 



The bicelle size varies depending on the initial con- 
ditions of the simulations. When a vesicle is used as 
the initial state instead of a random gas distribution, 
the ruptured vesicle opens up into a few large bicelles^i. 
The excess type A molecules form worm-like micelles or 
dissolve as isolated molecules. This scenario is differ- 
ent from the self-assembled bicelles in the present sim- 
ulations, where most of the type A molecules belong to 
bicelles for £ab < 1- Although the number of deter- 
gent molecules in bicelles can be changed via isolated 
molecules, the number of lipid molecules can be changed 
only by fusion and fission of micelles. Therefore, the 
growth of the bicelles becomes extremely slow after sta- 
ble bicelles, whose rims are completely covered by the 
detergents, are formed. The bicelle size is determined 
more by kinetics and depends on the initial states. Thus, 
the sample preparation method is important to control 
the bicelle size. 




FIG. 8: Sequential snapshots of vesicle rupturing into worm- 
like micelles at C* d = 1.2 and £ab = —5. (a) t = 7, 600to. (b) 
t = 8,200to. (c) t = 9,500r . (d) t = 30,000r . (a)-(c) Only 
the largest cluster is shown. The lower snapshot in (a) shows 
the cross-sectional view of the vesicle, (d) All molecules in 
the simulation are shown. 




FIG. 9: Time development of (a) the number N c \ us of 
molecules in the first and second largest clusters and (b) the 
aspericity a ap of the largest cluster for Cj^ d = 0.8 (•) and 1.2 
(x). The cluster size is normalized by the initial vesicle size, 
iVi n t = 2, 000. The same data are shown in Figs. [7] and [8] 
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FIG. 10: Detergent-adsorption-induced vesicle division and 
worm-like micelle formation, (a) Probability of structure 
changes at e* tt = 0.8 for C£ d = 0.8 (o), 1 (A), and 1.2 
(□). The solid lines represent the sum of vesicle division and 
micelle formation. The dashed lines represent worm-like mi- 
celle formation, (b) Phase diagram of vesicle structure for 
e£ tt = 0.6 (o), 0.8 (A), and 1 (□). The vesicle maintains its 
shape for eab values above the solid line. The vesicle is di- 
vided into two vesicles or transforms into worm-like micelles 
below the line. 



IV. SURFACTANT ABSORPTION INTO 
VESICLES 

In this section, we investigate the morphological 
changes of vesicles that are induced by detergent adsorp- 
tion. A lipid vesicle of type B molecules is placed in a 
solution of type A molecules [see Fig. EJa)]. The de- 
tergent adsorption induces vesicle division into two vesi- 
cles or rupture into worm-like micelles depending on the 
spontaneous curvature Cq (~ C^ d /2a) of the type A 
molecules and the attractive interaction between differ- 
ent molecules. 

Figures EHH] show typical examples of the vesicle divi- 
sion and rupture. The detergent adsorption first gener- 
ates an increase in the vesicle size [see Fig. Ufa)]. Subse- 
quently the vesicle transforms into a pear-shape, which 
leads to topological changes [sec Fig. EJb)]. The type A 
molecules are adsorbed onto the outer monolayer (leaflet) 
of the vesicle, and subsequently, they slowly move to the 
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FIG. 11: Time development of the mean number (A vcs ) of 
molecules in the vesicle, (a) sab = 0, —1, —2, —3, —4, and 
-5 at C*bd = 0.8. (b) C^d = 0, 0.2, 0.4, 0.6, 0.8, 1, and 1.2 
at eab = —4. Divided vesicles are not taken into account for 
the average in (b). 



inner monolayer. The resultant asymmetric distribution 
of the molecules induces an effective spontaneous curva- 
ture Cq S of the bilayer. The curvature Cg ff is caused by 
two mechanisms. One is the mismatch between the area 
difference of the two monolayers and the number differ- 
ence of the molecules in two monolayers as defined in the 
area difference elasticity (ADE) mode l 31 ' 60 ' 61 . The other 
is the spontaneous curvature of the bilayer, which is in- 
duced by the asymmetric distribution of type A molecules 
when Cq > 0. The outer monolayer contains more type 
A molecules than the inner monolayer. At Cq = 0, the 
curvature Cq is induced only by the area difference, and 
the contribution due to the molecular curvature increases 
with increasing Cq. 

For the vesicle division (see Fig. [7] and correspond- 
ing Movie£2), parts of the inner monolayer come into 
contact with each other on the inside of the pear neck 
and subsequently, the inner monolayer is separated by 
the outer (trans) monolayer (this saddle-point structure 
is called the trans- monolayer contact). The subsequent 
formation of a short cylindrical structure called stalk [see 
Fig- E£c)] leads to the division of the vesicle. The two 
resultant vesicles can have sizes equal to that of the orig- 
inal vesicle [see the red lines in Fig. Uta)]. This division 
pathway is the inverse version of membrane fusion in the 
modified stalk model^. While this pathway is dominant 
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FIG. 12: Vesicle initial growth rate fc a( j = (N vcs )ro/dt\t=o at 
C*bd = 0.4 (□), 0.8 (A), 1.2 (o) and eab = -4. The rate is 
normalized by the initial monomer concentration pi in (b). 



among the topological changes, other dynamics are also 
at play. At C£ d > 0.8, the large effective spontaneous 
curvature Cq S of the bilayer causes the opening of a pore 
at the neck of the pear-shaped vesicle. This in turn leads 
to the formation of worm-like micelles or to vesicle di- 
vision. The neck transforms into a cylindrical micelle 
structure and grows to form worm-like micelles. Other 
worm-like micelles grow from the edges of the flat bilayer 
(see Fig. [5] and corresponding Movie^). After the pore 
is opened, the surface area for absorption increases (the 
area of the inner monolayer and the worm- like arms), 
so that the absorption rate increases (see Fig. At 
C£ d < 0.6, no worm-like micelles are formed while some 
vesicles transform into worm-like micelles at C^ d > 0.8 
[see Fig. HPT a)]. In some simulations at C^ d > 0.4, 
the opened pores do not expand but lead to pinch-off of 
the neck [the cylindrical connection in the middle of Fig. 
E{a)] , and subsequently the two pores are closed resulting 
the formation of two vesicles (occasionally a vesicle and 
worm-like micelles). This pore opening is similar to the 
inverse version of the side-pore (or stalk-bending) fusion 
pathwa y 39 i 47 i 63 i 64 , although the present simulation does 
not exhibit the stalk intermediate stage during the divi- 
sion. At Cj^j < 0.2, the neck of the stalk structure often 
transforms to a bilayer disk and is detached as a disk-like 
micelle. 

The number of absorbed detergents increases with in- 
creasing attraction between the two types of molecules 



and with decreasing C^ d (see Fig. [TTj) . As eab decreases 
from to —3 for C^ d = 0.8, larger vesicles are formed. At 
£ab = —4, three of ten vesicles show topological changes 
and the other seven maintain a form of closed vesicles [see 
Fig. HOT a)]. The shape transition point is considered as 
the eab value at which half of the vesicles show either di- 
vision or rupture. The transition occurs at weaker levels 
of attraction for the middle range of C^ d ~ 0.6 as shown 
in Fig. HOT b). This is caused by the competition of two 
effects. With increasing C£ d , a lesser number of type A 
molecules are absorbed into the vesicle but the bilayer 
membrane becomes less stable for the same amount of 
absorption. We conclude that the detergents with mid- 
dle spontaneous curvature Cq ~ 0.3 (C^ d ~ 0.6) most 
efficiently induce vesicle division. 

The initial rate of vesicle growth k ad = {N ves )To/dt\ t= o 
is independent of eab when eab < — 1 [see Fig. [TTT a)]. 
since all the detergent molecules making contact with 
the vesicle are absorbed into the vesicle when the at- 
traction is sufficiently larger than the thermal energy. 
The growth rate k ad decreases with £^ tt as shown in Fig. 
[T2T a). This is caused by lower CMC values at larger 
values of e^ tt . A smaller amount of detergents remains 
as isolated molecules at larger values of e^ tt . When the 
growth rate fc a d is normalized by the concentration p\ of 
isolated molecules, it becomes almost constant, as shown 
in Fig. [T2Tbh If the concentration is higher than the 
CMC, the initial adsorption rate would not depend on 
Pa while the later evolution dynamics might be modified 
by pa- Here, we consider a fixed concentration of deter- 
gents pa- At higher values of pa, the divided vesicle may 
exhibit second and third divisions. 

Since a solvent-free model is employed, the vesicle can 
freely change its volume. This condition corresponds to 
high water permeability in comparison to the detergent 
adsorption rate. The theoretical study by Svetina and 
Bozi o 31 ' 32 suggests that higher permeability more eas- 
ily induces vesicle division. Our simulation results sug- 
gest that the addition of detergent molecules with slightly 
large head size efficiently induces vesicle division. 



V. SUMMARY 

We have studied structure formation dynamics in mix- 
tures of lipids and detergents. Bicelles are formed by 
self-assembly from a binary mixture of isolated molecules. 
The size of self-assembled bicelles increases with the in- 
creasing number fraction of lipids and repulsion between 
different molecules. When a vesicle is used as the an 
initial simulation state, a large bicelle can be formed, 
and the excess amount of the detergent molecules form 
worm-like micelles. In contrast, in self-assembly, most of 
detergent molecules are contributed to surround the rims 
of bicelles. Thus, the size of the bicelles are determined 
by the initial states and assembly kinetics. 

When a vesicle is placed in detergent solution, the de- 
tergent adsorption induces vesicle division or vesicle rup- 
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ture into worm-like micelles. The vesicle transforms into 
a pear shape and the contact of the inner monolayers 
(leaflets) on the inside of the pear neck induces vesi- 
cle division via the stalk intermediate in the modified 
stalk model. At large values of spontaneous curvature of 
the detergent monolayers, a pore is often opened in the 
pear neck, and this leads to vesicle division or worm-like 
micelle formation. A similar vesicle division has been 
experimentally observed in multi-component surfactant 
systems M i2!L3<L£L 

While our studies successfully reproduced two types of 
detergent-induced topological changes, other dynamics 
are observed in experiments such as the inside-out in- 
version of vesicles.— A pore-opening on the vesicle leads 
to the formation of a vesicle, whose outer monolayer is 
formed by the inner monolayer of the original vesicle. 
The adsorbed detergents form micelles with lipids in the 
outer monolayer and cause removal of the lipids from the 
outer monolayer, so that a negative spontaneous curva- 
ture of the bilayer is induced. This process has not been 
well understood and no simulations on this process have 



been reported. Since the size of spherical particles used 
for both molecules is the same in our simulations, they do 
not easily form micelles only from the outer monolayer. 
Thus, smaller detergent model molecules are likely nec- 
essary to simulate the vesicle inversion. The mechanisms 
to determine the assembly and lysis dynamics are not 
understood completely, and open problems remain for 
further studies. 
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